Capillary flow in an evaporating sessile drop 
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An analytical expression of velocity potential inside an evaporating sessile drop with pinned 
contact line is found. 

PACS numbers: 47.55.Dz 



I. INTRODUCTION 

The desiccated sessile drops attract the attention be- 
cause of different reasons and possible appUcations. ff a 
drop of pure Uquid dries on a smooth substrate, then the 
base of the drop shrinks but the contact angle remains 
fixed. Under evaporation, a drop of coUoidal suspension 
or solution having a strongly anchored three phase line 
keeps a spherical-cap shape with a constant base. The 
contact angle decreases with time and the height in each 
point of the profile decreases. To satisfy the anchoring 
condition, a flow of liquid has to occur inside the drop. 
An outward flow in a drying drop is produced when the 
contact line is pinned so that liquid that is removed by 
evaporation from the edge of the drop must be replen- 
ished by a flow of liquid from the interior P, 0, 0] . This 
flow is capable of transferring 100% of the solute to the 
contact line and thus accounts for the strong perimeter 
concentration of many stains. Ring formation in an evap- 
orating sessile drop is a hydrodynamic process in which 
solids dispersed in the drop are advected to the contact 
line. After all the liquid evaporates, a ring-shaped de- 
posit is left on the substrate that contains almost all the 
solute. Perhaps, everybody knows a dense, ring-like de- 
posit along the perimeter of a dried drop of coffee or tee, 
milk or juice on a table. 

Exploratory experiments Q using a variety of car- 
rier fluids, solutes, and substrates indicated that pref- 
erential deposition at the contact line is insensitive to 
a wide range of experimental conditions. Ringlike de- 
posits were observed whenever the surface was partially 
wet by the fluid irrespective of the chemical composition 
of the substrate. Different substrates were investigated 
(metal, polyethylene, roughened Teflon, freshly cleaved 
mica, ceramic, and silicon). Rings were found in big 
drops (15 cm) and in small drops (1 mm). They were 
found with aqueous and non-aqueous (acetone, methanol, 
toluene, and ethanol) solvents. They were found with 
solutes ranging in size from the molecular (sugar and 
dye molecules) to the colloidal (10 /Ltm polystyrene micro- 
spheres) and with solute volume fractions ranging from 
10~^ to 10""'^. Likewise, environmental conditions such 
as temperature, humidity, and pressure could be exten- 
sively varied without affecting the ring. Effects due to 



solute diffusion, gravity, electrostatic flelds, and surface 
tension forces are negligible in ring formation. 

The vertically averaged radial flow of the fluid in a des- 
iccated sessile drop was calculated in the work l]. The 
conservation of fluid was utilized. A constant evapora- 
tion rate all over drop free surface was assumed. 

For the case where the limiting rate is the diffusion of 
the liquid vapor, the evaporation of the drop rapidly at- 
tains a steady state so that the diffusion equation reduces 
to Laplaces equation Q . In this case the evaporation rate 
is larger near the contact line and the resulting ringlike 
deposit is more concentrated at the edge. 

The models Q and jsl deal with the vertically averaged 
radial flow of the fluid. The next Section describes how 
to obtain the space distribution of the fluid flow. 



II. MODEL AND RESULTS 

Likewise the works Q, H, |^, let us suppose, that 
the shape of the drop is a spherical cap. It means 
that the drop is small and a spherical cap shape in- 
duced by surface tension (Bond number smaller than one: 
Bo ~ g{p — pf)(Pa~^ <C 1, where g is gravitational ac- 
celeration, p is drop density, p/ is surrounding medium 
density , a is surface tension, d is drop diameter). The 
problem of a spherical cap on an impermeable substrate 
can be replaced by the problem of a lens. Let us con- 
sider an extremely simple situation when the cap is a 
semi-sphere, i.e. the contact angle is 90°. 

Under evaporation, the height in each point of the pro- 
file decreases. Nevertheless, under room condition the 
evaporation is a slow process. So, the desiccation time of 
a 15 mg sessile water drop is larger then 3500 s 3- This 
fact gives a possibility to consider a quasi-static process. 
We will consider a drop with a fixed free surface and 
the specific boundary conditions. Quasi-static approxi- 
mation was utilized in [3] to obtain a vapor rate near the 
free surface, too. 

Let us suppose that the evaporation rate is uniform all 
over drop free surface, i.e. we will consider only begin 
of the evaporation. The conservation of fluid determines 
the relationship between the velocity of the free surface, 
Un, the normal to the free surface flow of the fluid, w„, 
and the rate of mass loss per unit surface area per unit 
time from the drop by evaporation, J: 
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pVn + J = pUn- 
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Let us consider the potential flow, only. The poten- 
tial of the velocity field, iy9, is a solution of the Laplace's 
equation 



= 0. 



(1) 



Also, we will solve the Laplace's equation inside a spher- 
ical area using spherical coordinates . 

Measurements of the evolution of the drying drop 
height h with time t can be performed. Let us suppose, 
that the velocity of the drop height, uq, is known. Thus 
Un = uq cos 9, where is azimuthal angle. The volume of 
the drop is decreasing due to evaporation: 

(J — puo cos 6) dS — 0, 

where the integration goes all over the free surface. Thus 

_ uop 
~ 2 ■ 

The equation of solute conservation yields. 



r=R 



^.ol 2-|cos^|), (2) 



where R is the radius of the contact base. 

We looking for the boundary value problem Q 
inside a spherical area (r ^ R) as 



Taking into account Q equals ((SJ, if r = i?: 



^fc^i?'=-iPfe(cosi?) =;^6fePfc(cos^?). (6) 



fc=0 



R^ 



k=0 



Eq. El is always valid, if Ak — bkR/k. Hence, the 
coefficients of the series © are 

The capillary flow carries a fluid from the apex of the 
drop to the contact line (Fig. 
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where Pfc(cos'i?) are Legendre polynomials. 
Thus, 



FIG. 1; Capillary flow in a desiccated sessile drop 



fc=0 

Let us write the right hand side of the equation Q 
using Fourier-Legendre series expansion also known as a 
Generalized Fourier Series expansion: 

uo(^~\ cos^\ = hPkicosd). (5) 

Where, 

bk = ^^^"0 J Q - |cost9|^ Pfe(cosz9)sin^dz?. 



Thus, bo = 0, b2k+i = 0, 

b,, = uo{4k + l)22.(fc_i),(fc + i)r = 1> 2, . . . . 



III. DISCUSSION 

It is clear, more realistic assumption, that a drop has a 
spherical cap shape with an arbitrary contact angle, can 
not change the qualitative picture of the flow. In partic- 
ularly, the calculations of the vertically averaged velocity 
for the drops with the contact angles up to 90° Q did 
not show the qualitative changes of the flow. If a drop 
has a spherical cap shape with an arbitrary contact angle 
the Laplace's equation can be solved using toroidal co- 
ordinates. Although, the analytical results are too com- 
plicate in this case, and numerical calculations should be 
performed. Nonuniform evaporation rate can not change 
the results qualitative, but the line, where Vr{R) = 0, has 
to move to the base. 

Recently, the buckling instability was investigated dur- 
ing the drying of sessile drops of colloidal suspensions 
or of polymer solutions 0, S 0- Under solvent evap- 
oration, disperse particles or polymers accumulate near 
the vapor-drop interface. The outer layer of the drop 
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is more concentrated in the polymer and may display a 
gel or glassy transition and, licncc, may form a perme- 
able rigid gelled or glassy skin. This skin behaves like an 
elastic shell although it does not block the evaporation. 
This gelled or glassy skin will thus bend as the volume it 
encloses decreases, leading to large surface distortions. 



The processes of skin formation are not the scope of 
this brief communication. Nevertheless, knowledge of the 
velocity field inside an evaporating sessile drop is a neces- 
sary background to describe the colloidal particles trans- 
fer or solute diffusion, and skin formation. 
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